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Abstract 



We carry out a detailed numerical investigation of stochastic resonance 
in underdamped systems in the non-perturbative regime. We point out that 
an important distinction between stochastic resonance in overdamped and 
underdamped systems lies in the lack of dependence of the amplitude of the 
noise-averaged trajectory on the noise strength, in the latter case. We provide 
qualitative explanations for the observed behavior and show that signatures 
such as the initial decay and long-time oscillatory behaviour of the temporal 
correlation function and peaks in the noise and phase averaged power spectral 
density, clearly indicate the manifestation of resonant behaviour in noisy, 
underdamped bistable systems in the weak to moderate noise regime. 
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I. INTRODUCTION 

Noise usually has a disruptive role in na- 
ture. However, there is a phenomenon [1,2] 
in which noise helps in making the system 
behave in a more coherent manner. This 
intriguing and rather counterintuitive phe- 
nomenon called Stochastic Resonance(SR) [3] 
was first proposed as an explanantion of the 
observed periodicity in the ice ages on earth 
[1,2], but has since been observed in a large 
variety of physical, chemical and biological 
systems [4]. Bistable systems subject to noise 
and periodic modulation constitute the most 
common class of systems which exhibit SR. 

The phenomenon of SR can be thought 
of as resulting from the synchronization of 
a stochastic time-scale (determined by the 
thermal transition rate over the barrier) and 
a deteministic time scale (determined by the 
time-scale of periodic modulation). When 
the noise strength is appropriately tuned to 
ensure the synchronization of the these two 
time scales, the system exhibits stochastic 
resonance. SR manifests itself through a va- 
riety of signatures such as peak in the am- 
plitude of the noise-averaged trajectory ver- 
sus noise strength, peak in the distribution 
of peak heights in the residence time dis- 
tribution (RTD) [5-7], peaks in the power 
spectral density (or Signal-to- Noise Ratio) [8] 
at the modulating frequency and late-time 
oscillatory behaviour of the temporal auto- 
correlation function [5]. SR has been exten- 
sively investigated in many overdamped sys- 
tems, but it has'nt been explored in great 
detail in underdamped systems, (references 
listed in [5,9] are some exceptions). Apart 
from the intrinsic interest in exploring in de- 
tail the manifestations of SR in the under- 
damped regime, there are other motivations 
for focussing on this regime. In particular, 
relativistic systems are governed by under- 
damped equations of motion and investigat- 
ing SR in single particle underdamped sys- 



tems may provide valuable insights into the 
manifestations of SR in relativistic field the- 
ories. 

We will focus on numerically investi- 
gating the phenomenon of SR in under- 
damped bistable systems subject to white 
noise and periodic modulation. The system 
we are dealing with can be described by the 
Langevin equation 

d^cc dec 

-jw + V-tt - x + x 3 = A cos(ujt + 9)+ f (t) 
at z at 

{t(t)W)) = ZTr,5(t-lf) (1) 

where T is the temperature of the back- 
ground heat bath, r\ is the dissipation co- 
efficient, A Q and uj are the amplitude and 
frequency of periodic modulation repectively. 
£ is a Gaussian distributed (delta corre- 
lated) white noise term which satisfies the 
fluctuation-dissipation theorem. For simplic- 
ity, we have made use of scaled, dimensionless 
quatities in the above equation. The poten- 
tial governing particle dynamics is of stan- 
dard double well type, modulated by a peri- 
odic signal with small amplitude. 

1 1 

V(x, t) = — x 2 + -x A - A x cos(ut + 6) (2) 

In the absence of periodic modulation, the 
doubly degenerate minima of the double well 
potential are located at x m = ±1 and the bar- 
rier height is given by AVo = 0.25 The modu- 
lation amplitude Aq is small enough to ensure 
the bistability of the potential at all times. 
Moreover, the noise strength, D = Trj (for 
the underdamped case) is much smaller than 
the barrier height in the absence of modula- 
tion i.e. D « AVo. We will also be working 
in the moderately weak noise regime D ^$ A . 
There has been several analytical studies of 
SR in overdamped systems [3], but those in- 
vestigations were based on perturbation the- 
ory and therefore restricted to the perturba- 
tive regime Aq << D or the regime of va- 
lidity of the adibatic approximation A « 
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D « AVo- An analytical study of un- 
der damped, periodically modulated systems 
was carried out by Jung and Hanggi [10], 
based on the solution of the 2-dimensional 
and 3-dimensional Fokker-Planck equations 
for probability distribution. They obtained 
the asymptotic probability distribution and 
showed that for uniformly distributed initial 
phases, the power spectrum exhibits delta- 
functions spikes at multiples of the driving 
frequency. Previous experimental work on 
SR in underdamped systems [5] have been 
restricted to the strong noise regime char- 
acterized by Trj >> AV and focussed pri- 
marily on signatures like the signal-to-noise 
(SNR) ratio and residence time distribution. 
More recently, Alfonsi et.al. [9] have shown 
that that an intra-well SR can coexsist with 
conventional inter-well SR in underdamped 
systems and is manifest as twin peaks in the 
plot of power spectral density versus noise 
strength. Their paper further extended the 
work of Stocks et.al. dealing with SR in 
monostable systems [11] using methods of 
Linear Response Theory. Evidence of SR in 
overdamped bistable systems without explicit 
symmetry breaking and subject to periodi- 
cally oscillating noise strength have been dis- 
cussed in [12,13]. 

Our main aim in this paper is to numeri- 
cally explore the various manifestations of SR 
in underdamped systems in the moderately 
weak noise and moderate forcing regime. 
This is distinct from the strong noise, weak 
forcing regime [5] and strong forcing regime 
[9] explored in earlier works. Although the 
effect of weak noise was considered in the 
work of Alfonsi et.al. [9], their analysis was 
focussed on the behaviour of the power spec- 
trum and in distinguishing between intrawell 
and interwell SR in the high frequency forc- 
ing regime. We point out an unusual be- 
haviour of the noise- averaged trajectory, in 
under-damped systems in the mo derate- weak 
noise and weak forcing regime, which is quite 



distinct from its behaviour in over-damped 
systems where it is widely used charac- 
teristic signature of SR. We find that there is 
no increase in the amplitude of noise averaged 
trajectory with increasing noise strength, un- 
like SR in overdamped systems. This im- 
portant distinction in the dynamics of over- 
damped and underdamped systems was not 
recognized earlier. It can be attributed to 
the delicate balancing act between noise and 
dissipation in underdamped systems, which 
also reduces the window size in parameter 
space for observing SR-like behaviour. Our 
results supplement the analysis of SR in un- 
derdamped systems carried out earlier [5,9]. 

The paper is organized as follows. In Sec. 
II, we carry out a detailed numerical study 
of SR in underdamped systems in the weak 
noise, moderate forcing regime. The effect of 
changing system parameters on the various 
distinct signatures of SR are described in de- 
tail. We conclude with a brief summary and 
discussion of our results in section III. 

II. STOCHASTIC RESONANCE IN 

UNDERDAMPED SYSTEMS : 
MODERATE AND WEAK NOISE 
REGIME 

The key idea behind SR is that the syn- 
chronization of the stochastic and determin- 
istic time scales of the system can result 
in a more coherent response of the system 
even on increasing the noise strength, for a 
small range of values of the modulating fre- 
quency. Even though SR was first inves- 
tigated for overdamped systems, this intu- 
itive notion of SR was also found to be valid 
for underdamped systems [5]. The interest- 
ing area of parameter space is such that nei- 
ther the periodic modulation of the potential, 
nor the thermal noise activation is individ- 
ually capable of ensuring periodicity in the 
inter-well hopping of the particle. For large 
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noise strengths i.e. D » AVq the particle 
has a large enough thermal energy to easily 
surmount the potential barrier and therefore 
does not feel the presence of the barrier. In 
the absence of coupling to a heat bath, the 
small amplitude of modulation, Aq < AVo 
and forcing frequency tv << r^ 1 (where r r 
is the intra-well relaxation time scale) can- 
not provide sufficient energy to the parti- 
cle to surmount the barrier. However, both 
these effects acting in conjunction can, under 
certain circumstances, lead to periodic be- 
haviour in the noise-averaged particle trajec- 
tory. In general, the thermal transition rate 
(Kramer's rate) over the barrier r^T 1 = r k , 
depends on both the noise strength and dis- 
sipation coefficient. The time scale of pe- 
riodic modulation (r = -j) is determined 
by the forcing frequency uj. Due to the odd 
power of the coupling between the particle 
position and the periodic forcing, the symme- 
try of the potential is explicitly broken. The 
particle in either well encounters the small- 
est (largest) barrier twice during a single pe- 
riod. If (t^) are the corresponding 
thermal transition time-scales, the noise av- 
eraged particle trajectory becomes periodic 
when t\ ■* < tq < T k + \ Inter- well hopping 
is therefore synchronized when the following 
(approximate) condition between the deter- 
ministic and stochastic time scales is satisfied 
r k ~ 2> =^ Usr ^ jl # p ox an( j l u [14] have 

shown that this matching condition which 
yields the corresponding resonant value of the 
noise strength (D sr = (Tr)) sr ), is an approxi- 
mate one and does not become exact in any 
limit. 

At this stage it is necessary to point out 
the crucial difference between our results and 
those obtained via analog simulations [5] in 
the work of Gammaitoni et.al. The main dif- 
ference lies in the distinctly different regime 
of study undertaken in this paper. Even 
though the work of Gammaitoni et.al. [5] 



dealt with investigating SR in underdamped 
systems, their choice of parameters restricted 
their study to the strong noise regime charac- 
terized by the condition Tr\ > AVq. However, 
we decided to explore the weak-to-moderate 
noise regime characterized by the condition 
Trj << AVq and Trj < Aq. Moreover, the va- 
lidity of their theoretical analysis depended 
on the applicabiliy of perturbation theory, 
which amounts to the condition A « Trj 
implying that the influence of modulation is 
small compared to thermal activation. This 
lead to the temperature dependence of the 
amplitude of the noise-averaged trajectory. 
However, our regime of study is not amenable 
to a perturbative analysis since in our case 
Trj < Aq. Moreover, Aq < AVo further im- 
plying a substantial (non-perturbative) mod- 
ulation of potential by the driving force, even 
though the modulation is insufficient to de- 
stroy the bistability of the potential at any 
time and is therefore distinct from regime 
where resonant trapping [15] is observed. We 
find that this leads to an important difference 
in the behaviour of the noise averaged trajec- 
tory characterized by the lack of dependence 
of the amplitude with noise strength. This is 
distinct from the manifestation of SR in the 
overdamped case and in the underdamped 
case studied earlier [5]. 

The particle dynamics in the moderate- 
strong damping regime is dominated by fluc- 
tuations and the Kramer's rate in this regime 
is given by 

However, this expression is derived under 
the assumption that the intrawell relaxation 
time-scale of a particle is much smaller than 
the thermal transition time scale and the par- 
ticles are able to thermalize before escaping 
over the barrier. As the dissipation coefficient 
decreases, this assumption breaks down as 
the motion of the particles starts being domi- 
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nated by energy diffusion. Appropriate mod- 
ification to the rate formula in the very weak 
dissipation regime was obtained by Kramers 
[16,17], Carmeli et.al. [18] and by Buttiker 
et.al. [19], in which the transition rate was 
found to be proportional to i] 



(4) 



where I(E ) is the action at the top of the 
barrier. [Carmeli et.al. [18] and Buttiker 
et.al. [19] also obtained appropriate correc- 
tions to the above formula by relaxing the 
approximation that the distribution function 
drops to zero at the top of the barrier, a cru- 
cial assumption in the derivation of Eq.(4).] 
The approximate resonant frequency uo sr ~ 
^- for underdamped systems has to be esti- 
mated with this modified rate formula in the 
very weak dissipation regime. In our work, 
uj sr was used only to identify the regime for 
resonant behaviour in underdamped systems. 

We numerically solved Eq.(l)using a sec- 
ond order stochastic Runge-Kutta algorithm. 
The noise averaged trajectory was obtained 
by averaging over one hundred noise realiza- 
tions. The upper plot of Fig.l shows the noise 
averaged trajectory for two different values 
of bath temperature, obtained by solving the 
overdamped Langevin equation given by 



^ - x + x 3 = A cos{ut) + £(t) 
at 

mm) = 2TS(t-t>) 



(5) 



A clear dependence on the amplitude of the 
oscillation on temperature is observed. The 
maximum of the amplitude corresponds to 
the temperature that satisfies the resonant 
condition u sr ~ — . When the tempera- 

T k 

ture increases beyond 0.08, the amplitude de- 
creases and eventually the dynamics becomes 
noise dominated resulting in loss of period- 
icity. In contrast, the lower plot of Fig.l 
shows the noise averaged trajectories for un- 
derdamped case obtained by solving Eq.(l). 



The trajectories for the two different temper- 
ature values almost overlap indicating a lack 
of dependence on the temperature of the heat 
bath. This important aspect of SR in under- 
damped system in the weak-moderate noise 
regime was not realized earlier. It can be at- 
tributed to the fact that the noise strength 
(Trj) in underdamped systems depend not 
only on the temperature of the heat bath, but 
also on the dissipation coefficient r\. Since r\ 
is very small for underdamped systems, an 
increase in temperature is offset by the small- 
ness of rj which therefore does not lead to any 
substantial increase in noise strength. 



A. Residence Time Distribution 

An useful characterization of SR [5] is ob- 
tained by considering the discrete stochastic 
process associated with the barrier crossing 
time. The residence time then corresponds 
to the time spent by the particle in a well be- 
tween successive transitions over the barrier 
and is obtained by taking the difference be- 
tween successive barrier crossing times. The 
residence time was obtained numerically by 
noting the time required by the particle start- 
ing from either well to cross a point specified 
on the other side of the barrier for each tra- 
jectory. The average of these times over all 
realizations yielded the mean residence time. 
The choice of the crossing point, although ar- 
bitrary, was made to ensure that the particle 
had indeed made a transition from one well 
to another. This rules out trajectories where 
the particle just crosses the barrier, lingers 
around the top of the barrier, before falling 
back into the well it started from. This phe- 
nomenon is particularly relevant in under- 
damped systems and would lead to under- 
estimating the residence time. We empha- 
size that the two-state approximation was not 
been made in estimating the residence time. 
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Although it is often customary to make such 
an approxiation in analyzing hopping dynam- 
ics in overdamped systems, it would lead to 
wrong estimates of the residence time for rea- 
sons mentioned above. Fig. 2 shows the resi- 
dence time distribution for 3 different values 
of the bath temperature. A pattern of peaks 
in the RTD with exponentially decaying peak 
heights is observed. (The peak height at 
tjr = 0.5 have been normalized to unity 
in Figs. 2,4, 6.) The likelihood of the particle 
hopping over the barrier is largest when the 
barrier height to be surmounted is the least. 
Due to the modulation of the potential, the 
particle encounters the smallest barrier twice 
during every modulation period. When per- 
fect synchronization is achieved i.e. the res- 
onant condition is satisfied, the particle hops 
over the barrier twice every modulation pe- 
riod and a single large peak in the RTD is ob- 
served at t — 0.5t . When T > T sr (where T sr 
is the bath temperature which corresponds to 
the resonant frequency u sr for fixed 77) , noisy 
dynamics starts dominating and synchroniza- 
tion of inter- well hopping is lost. The inter- 
well transitions occur for t < 0.5r , which ac- 
counts for a sharp peak in RTD for t < 0.5ro. 
For T < T sr , perfect synchronization is lost 
because the thermal transition time scale(Tfc) 
is larger than 0.5r . If the particle in un- 
able to traverse the barrier when the barrier 
height is minimum, it has to wait for one full 
period before it again faces the smallest bar- 
rier. This accounts for the multiple peaks in 
the RTD at odd multiples of 0.5r . 

We would like to emphasize that the role 
of dissipation on SR can only be addressed by 
studying UD systems. This is easily under- 
stood since the thermal transition time scale 
for inter-well hopping depends on the 77, al- 
beit weakly. However, the signature of SR 
associated with the periodicity of the noise- 
averaged trajectory is not very sensitive to 
changes in r\ and uj. This can be seen from 
the plots of the noise averaged trajectory for 



different values of 77 and uj as shown in Fig. 3 
and Fig.5. The periodicity in the noise av- 
eraged trajectory is observed even when 77 
is increased by about two orders of magni- 
tude. A more sensitive indicator of reso- 
nant behaviour is the residence time distri- 
bution. Nevertheless, the dependence on 77 
is rather weak as is evident from the plots 
for the RTD shown in Fig.4. Increasing 77 
increases the noise strength thereby increas- 
ing fluctuations, but that increase is balanced 
out by the damping of fluctuations resulting 
from a larger rj. Because of this delicate bal- 
ance between fluctuation and dissipation ob- 
served in underdamped systems, the depen- 
dence of RTD on 77 is weaker than its de- 
pendence on the bath themperature. This 
is evident from Fig.4 which shows plots of 
the RTD for three different values of 77. For 
very small dissipation, r k becomes larger than 
0.5to and this leads to loss of synchronization 
of inter-well transitions resulting in multiple 
peaks in RTD at odd multiples of 0.5ro. How- 
ever, the weak dependence on 77 does not af- 
fect the shape of the RTD on increasing 77 
by an order of magnitude (compared to the 
benchmark value of 77 = 0.03) as is evident 
from the last two plots of Fig.4. 

Fig.5. shows the the noise averaged tra- 
jectory for three different values of the mod- 
ulating frequency. The periodic behaviour is 
observed even for values of uj substantially 
different from the value which satisfies the 
resonant condition. However, as uj is in- 
creased, synchronization is lost and the dy- 
namics becomes noise dominated as is evi- 
dent from the last plot of Fig.5. However 
the sensitivity to change in uj is more pro- 
nounced in the plots for the RTD shown in 
Fig. 6. When perfect synchronization is at- 
tained, inter-well hopping over the barrier oc- 
curs when the the barrier height is minimum. 
This amounts for a single sharp peak in the 
RTD at t — 0.5r corresponding to the reso- 
nant frequency uj sr = 0.03, as seen in the first 
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plot of Fig. 6. For larger (or smaller) mod- 
ulation frequencies, synchronizations is lost 
and this is manifest through appearance of 
multiple peaks (with exponentially decreas- 
ing peak heights) in the RTD at odd multi- 
ples of the half-period of modulation. 

However, a proper signature of SR based 
on RTD is characterized by a peak in the dis- 
tribution of peak heights (at half the driv- 
ing period) versus noise [6]. We also find a 
peak around T sr , in the distribution of peak 
heights plotted against temperature. How- 
ever, no such peak is observed when the dis- 
tribution of peak heights are plotted against 
the dissipation coefficient, providing another 
indication that the hopping dynamics is less 
sensitive to dissipation coefficeint in the weak 
damping regime. This characterization of 
SR has lead to some controversy regarding 
the identification of SR bonafide reso- 
nance [7,20-22] in the overdamped case. Choi 
et.al. [7] have argued that for large A /D val- 
ues, peaks in the RTD distribution were in- 
sufficient to imply resonant behaviour in a 
practically viable frequency range. Our aim 
in carrying out the RTD analysis for vary- 
ing temperature, dissipation and driving fre- 
quency was to indicate that the behaviour 
of the RTD in underdamped systems follows 
the same pattern seen in overdamped systems 
and can be understood in terms of synchro- 
nization of hopping over the barrier. How- 
ever, since we are working in the large A /D 
regime, RTD alone cannot be used to in- 
fer about the manifestation of SR in under- 
damped systems. In order to unambiguously 
show that SR is observed in underdamped 
systems in the weak and moderately weak 
noise regime (i.e.A /(Tr)) >> 1) investigated 
in this paper, we present below two other sig- 
natures of SR. 



B. Auto-correlation Function and Power 
Spectrum 

Evidence for SR can also be obtained 
by studying the behaviour of the auto- 
correlation function (ACF) [5] , defined by the 
equation 

C(t) = (x(t)x(0)) (6) 

where the (...) implies averaging over noise 
realizations. The first plot in Fig. 7 shows the 
early-time exponential decay of the ACF for 
four different values of the bath temperature. 
Larger bath temperatures lead to a faster de- 
cay in the ACF, C(t), since large fluctua- 
tions destroy correlations at large temporal 
separations. The initial decay of the ACF is 
followed by oscillatory behaviour with a fre- 
quency equal to that of the the modulation 
frequency. Such a behaviour is symptomatic 
of SR in driven, diffusive systems subject 
to thermal noise. The lower graph of Fig. 7 
shows a distribution of the temporal ACF at 
late times which indicates a clear oscillatory 
dynamics. The gray solid line is a sinosoidal 
fit to the numerical ACF data and has a fre- 
quency equal to the modulation frequency of 
the system. Note that our choice of parame- 
ters does not allow for a perturbative analysis 
which was carried out in [5] to obtain anayti- 
cal expressions for the amplitude of the noise 
averaged trajectory and the auto-correlation 
function. Such an approximation is only valid 
for small perturbations A Q x m « T. 1 How- 
ever, in our case, T > AVo- Furthermore, our 
regime of study lies beyond the perturbative 
regime, since A x m ^ T. 

Another common signature of SR can be 
obtained by looking at the noise and phase 



In the notation of [5], D corresponds to the temperature (T) of the heat bath. 
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averaged power spectral density S(u) ob- 
tained by taking the fourier transform of the 
ACF C(t). 

/oo 
exp"^* << x(t + r)x(t) » dr 
-oo 

(7) 

where the inner brackets indicate avergaing 
over noise realizations and the outer brack- 
ets denote avergaing over the initial random 
phase. We averaged over 50 different noise 
realizations and 20 different values of ran- 
dom initial phase to first obtain the noise and 
phase averaged temporal ACF for 64 different 
values of temporal separation. The fourier 
transform of the noise and phase averaged 
ACF then yielded the noise and phase aver- 
aged power spectral density. Fig. 8 shows the 
power spectral density function of the 
frequency, obtained for four different values 
of the bath temperature, all other parame- 
ters remaining fixed. The plots clearly show a 
sharp peak at the the modulation frequency, 
and indicates that the peak strength at the 
driving frequency goes through a maximum 
as a function of the noise strength. Although, 
the exact location of the maxima is difficult 
to ascertain precisely, it is nevertheless clear 
that the maxima is close to T = 0.2, which 
is the temperature at which the synchroniza- 
tion condition is satisfied for the given choice 
of parameters. These additional signatures 
clearly establish that the observed synchro- 
nization, in the weak and moderately weak 
noise regime of underdamped bistable sys- 
tems, is indeed a manifestation of stochastic 
resonance. 

III. CONCLUSIONS 

To summarize, we have carried out a 
numerical study of SR in underdamped, 
bistable systems in the regime characterized 
by the conditions Ti] « AV and Tr] « 



Aox m . This is distinct from the regime in- 
vestigated earlier [5] and is not amenable to 
an analysis based on treating the modulation 
of the potential as a small perturbation. A 
comparison between overdamped and under- 
damped dynamics was carried out. No sig- 
nificant dependence of the amplitude of the 
noise averaged-trajectory, on the tempera- 
ture of the heat bath, was found in the under- 
damped scenario. This is in marked contrast 
to the the overdamped case, where the am- 
plitude depends sensitively on the bath tem- 
perature. Synchronization of hopping with 
periodic modulation of the potential well was 
observed in the behavior of the RTD, as in the 
overdamped case. The RTD was found to be 
more sensitive to changes in bath tempera- 
ture and only weakly dependent on changes 
in the dissipation coefficient. The auto- 
correlation function and the power spectral 
density were also obtained for various val- 
ues of the bath temperature to unambigu- 
ously confirm that the underdamped system 
does indeed exhibit resonant behaviour in the 
weak noise and weak forcing regime. The 
power spectrum was characterized by a sharp 
peak at a value close to that of the modulat- 
ing frequency, with the largest peak height 
corresponding to the temperature at which 
the resonant condition was satisfied. These 
signatures provide evidence that SR in un- 
derdamped systems depends on the interplay 
between the bath temperature and the dis- 
sipation coefficient and leads to a rather dis- 
tinct pattern of dynamics, not observed in the 
previously investigated regimes of parameter 
space. 

An interesting extension of this work 
would involve looking for evidence of SR in 
relativistic field theoretic models. The study 
of stochastic resonance for spatially extended 
systems has been carried out for Ginzburg- 
Landau type field theories [23], but has been 
restricted to the over-damped regime .There 
it was found that an appropriate choice of the 
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frequency of the periodic driving obtained by 
matching the thermal activation time- scale 
to half the period of the modulating back- 
ground, can result in periodically synchro- 
nized behavior of the mean field about = 
(see Fig. 2 of Ref. [23]). It would be inter- 
esting to explore the possibility of observ- 
ing SR in underdamped, field-theoretic mod- 
els exhibiting spontaneous symmetry break- 
ing of a Z 2 symmetry since it would provide 
the first evidence of emergence of coherent 
behaviour in noisy, underdamped field the- 
ory. 
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OD Dynamics 
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UD Dynamics 




FIG. 1. Plots showing the noise averaged trajectory with two different temperatures for the 
overdamped case (top) and underdamped case (bottom). The bath temperatures for the two tra- 
jectories shown correspond to T = 0.04,0.08 (top) and T = 0.08,0.2 (bottom). 
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FIG. 2. Residence Time Distribution for three different values of bath temperature. Other 
fixed parameter values are r\ = 0.03, Aq = 0.1, u = 0.03. The temperature at which the resonance 
condition is satisfied for fixed oj and 77 is T = T sr = 0.2 
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FIG. 3. Noise averaged trajectory for three different values of rj. Other fixed parameter values 
are T = 0.2, Aq = 0.1, to = 0.03. The dissipation at which the resonance condition is satisfied for 
the given w and T value is rj = i] sr = 0.03 
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FIG. 4. Residence Time Distribution for three different values of r/. Other fixed parameter 
values are T = 0.2, A = 0.1, w = 0.03. 
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-0.4 

FIG. 5. Noise averaged trajectory for three different values of the modulating frequency. Other 
fixed parameter values are T = 0.2, Aq = 0.1, rj = 0.03. The resonant frequency for the given r\ and 
T values correspond to uj = uj sr = 0.03 
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FIG. 6. Residence Time Distribution for three different values of the modulating frequency. 
Other fixed parameter values are T = 0.2, Aq = 0.1, rj = 0.03. The resonant frequency for the given 
r] and T values correspond to w = uj sr = 0.03 
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FIG. 7. The upper plot shows the evolution of the temporal auto-correlation function (ACF) 
at early times and clearly indicates the exponential decay of the ACF. The four curves correspond 
to T = 0.2, 0.4, 0.6, 0.8 respectively from top to bottom. The lower plot shows the oscillatory 
behaviour of the ACF at late times for T=0.4. The solid line is a sinosoidal fit to the numerical data 
with an oscillating frequency to = 0.03. Other parameter values are r\ = 0.03, a; = 0.03, Aq = 0.1. 
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FIG. 8. Plots of the power spectrum showing peaks at the resonant frequency uj st = 0.03. The 
largest peak is seen for T = T sr = 0.2 which corresponds to the temperature at which the resonance 
condition is satisfied for uj = 0.03, rj = 0.03, Aq = 0.1. 
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